Universal scaling in nonequilibrium transport through an Anderson impurity 
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Using non- equilibrium renormalized perturbation theory, we calculate the conductance G as a 
function of temperature T and bias voltage V for an Anderson model, suitable for describing trans- 
port properties through a quantum dot. For renormalized parameters that correspond to the extreme 
Kondo limit, we do not find a simple scaling formula beyond a quadratic dependence in T and V. 
However, if valence fluctuations are allowed, we find excellent agreement with recent experiments. 
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Universality is one of the most beautiful and useful 
concepts in physics. In general, the physical proper- 
ties of a system depend on a certain number of parame- 
ters which change for different experimental realizations. 
However, in favorable cases, physical observables are de- 
scribed by the same universal function, once the different 
physical magnitudes are scaled appropriately. For exam- 
ple in the case of the temperature dependence of the con- 
ductance through one quantum dot G(T) in the limit of 
zero bias voltage V, once a characteristic energy scale Tk 
(the Kondo temperature) is identified, the conductance 
of different systems is very well described by the same 
universal function G(T /Tk), even if the systems have 
very different Tk*&> Scaling and universality are con- 
cepts which are quite naturally connected to renormal- 
ization group treatments of the Anderson model in the 
Kondo regime (Coulomb repulsion U much larger than 
the resonant level width A), and in fact numerical renor- 
malization group (NRG) calculations reproduce the scal- 
ing mentioned above and in other physical properties. 

Theoretically, the situation is much more difficult in 
the nonequilibrium situation which arises for a finite 
bias voltage between the leads connected to the quan- 
tum dot in the experiment. Only recently, extensions to 
the nonequilibrium case of essentially exact techniques 
such as NRG 5 and exact Bethe ansatz^ were proposed, 
while approximations used at equilibrium have shortcom- 
ings when extended to the nonequilibrium c&seZ Never- 
theless, using a Fermi liquid approach, based on pertur- 
bation theory (PT) in U/A, and Ward identities, Oguri 
has determined exactly the scaling for T and eV small 
compared to Tk for the Anderson model 8 
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where Go = G(0, 0) and the values of ct and a are dis- 
cussed below. 

Recent experiments in GaAs quantum dots for differ- 
ent situations in the nonequilibrium regime,— for low T 
and V have found that G(T, V) is well described by a 
universal scaling function that extends Eq. (|TJ) to higher 
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Here c T ^ 5.488 is fixed by Eqs. ((TJ and ©, a = 0.10 ± 
0.015, 7 = 0.5 ± 0.1 and Ge(T) is an empirical curve 
obtained from a fit to NRG results: 
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with s = 0.21 for an impurity with total spin S = 1/2. 
From these equations, one can see that a is the ratio of 
the term of order [eV/(Jfc7jc)] 2 with respect to that of 
order (T/Tk) in the decrease in the conductance, while 
7 represents the effect of terms [eV/(kT K )} 2 (T/T K ) 2n 
with integer n > 1. 

From an exactly solvable anisotropic Kondo model^ 
one extracts a = 3/tt 2 ~ 0.30 and 7 = 2(ttT k /T a ) 2 / 'c T ~ 
3.60(Tft7T a ) 2 , where T a is an energy scale of the order 
of Tk- To our knowledge, no other precise information 
on 7 exists. 

The purpose of this work is to test the observed scaling 
relation and calculate 7 in the impurity Anderson model, 
using renormalized PT (RPT)J^ The basic idea of RPT 
is to reorganize the PT in terms of fully dressed quasi- 
particles in a Fermi liquid picture. The main advantage 
is that even in the strong coupling (SC) limit U —> 00, 
for which ordinary PT in u = U/(ttA) becomes invalid, 
the corresponding ratio between renormalized parame- 
ters (denoted by a tilde) becomes u = U/(irA) — > 1, 
being u < 1 for finite U. 1 ® For nontrivial cases, already 
free quasiparticles (taking [7 = 0, which is similar to slave 
bosons in a mean field approximation 12 ) reproduce very 
well the low-frequency part of the equilibrium spectral 
density at the quantum dot. An example is a case in 
which the Kondo peak is split in twoJ^ U (proportional 
to the full vertex) represents the "residual" interaction 
between quasiparticles. Calculating the renormalized re- 
tarded self-energy E r within nonequilibrium RPT up to 
second order in u, T and V leads to the exact result Eq. 
(JTJjSiii Here we calculate numerically S r up to order u 2 , 
for finite kT and eV but smaller or of the order of the 
Kondo energy kTjc- 
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Ordinary PT up to second order in U supplemented 
by an interpolative perturbative approach (IPA) ; 13 i 14 
(which corrects the second-order result in order to re- 
produce exactly the atomic limit U / A — > +00) has been 
shown to describe well the conductance through a quan- 
tum dot for U < 8A.— The results agree with those ob- 
tained using the finite temperature density matrix renor- 
malization group method^ Comparison of the spin de- 
pendent IP A 17 i 18 with exact diagonalization in finite sys- 
tems shows very good agreement for U = 6.25A»i£ The 
extension of PT in U 2 to the noncquilibrium case has 
been considered by Hershfield et al*^ They found that 
for finite V, the current is conserved only in the electron- 
hole symmetric Anderson model (SAM). Different self- 
consistent approaches were proposed to overcome this 
shortcoming, by a suitable election of the unperturbed 
Hamiltonian.— ^ While these approaches work well in ab- 
sence of a magnetic field B, numerical difficulties per- 
sist for small non-vanishing B and Here we will take 
B = and parameters corresponding to the SAM for 
our numerical integrations. In this case the current is 
conserved for each spin without the need to solve self- 
consistent equations^ 

We use the spin 1/2 Anderson model to describe a 
quantum dot interacting with two conducting leads, one 
at the left and one at the right, with chemical potentials 
/jLl and /j,r respectively, with jjll — fiR = eV. We define 
the zero of energy by [II — sVAr/A, where A = A^ + 
A R and A„ = irj^k \Vku\ 2 S(uj - e a eff ) (neglecting here 
the smalt dependence on to). The Hamiltonian is split 
into a noninteracting part Hq and a perturbation H' as 

H = H + H', 

H = ^2e k „ c kva c kva + ^2 e%ff n da 

kver 

H' = (Ed - elff) n da + U n d -in dll 
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where v = L,R refers to the left and right leads. In 
general e a e ^ should be determined selfconsistently, but 
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conductance G = dl/dV from numerical differentiation 
of the current /, which can be written as 20 
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where f u (u) = f(u - fi u ), /(w) = l/{e"' kT + 1), 
A = AAlAr/ A 2 indicates the degree of asymmetry 
of the hybridization of the dot with both leads, and 
p{uj) = — ImG da (uj)/iT where G dcr (uj) is the retarded 
Green's function of the electrons at the dot for spin c, 
which can be written asi 
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Within RPT, the low frequency part of G d(7 (uj) can be 
approximated as^ 
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where z = [1 - dT?Jdw]-\ e° ff = z[e° ff + E£(0)], 

A = zA, and £™ m M = z^ a cm (ui) , where the remain- 
der retarded self-energy is defined as 

£? m M - - K(P) ~ udrrjdu. (8) 

In Eqs. and ©, E£(0) and dT^Jdw are evaluated 
at uj = T = V = 0. A comparison between G d(T (uj) (cal- 
culated within PT) and G da (ui) with E£(w) = 0, for a 
case with non-trivial frequency dependent A(w) is pro- 
vided in Ref. |12| . showing a very good agreement for low 
For large values of U/A, ordinary PT in U is not re- 
liable and z is in principle not known, although it can be 
obtained from exact Bethe ansatz calculations. However, 
replacing A by A/z in Eq. ([5]), z cancels and the current 
is expressed in terms of renormalized parameters £gjj, 

A, and E™ m (w). In the SC limit U/A -> 00, Hewson 10 
has shown that the ratio of renormalized parameters is 
u = U/(ttA) = 1 and furthermore, defining the Kondo 
temperature Tjf by the linear term in the specific heat 
in this limit 7c = 7r 2 fc/(6T£), one obtains A = 4kT£/ir. 
In the SAM at B = 0, which we shall use, A L = A R , 
E d = —U/2, and this implies £°t* — 0; 12 Experimen- 
tally, Al ~ A R and the scaling results do not depend on 
the gate voltage, which controls E d . Furthermore, the 
value of E d is irrelevant in the SC limit. This justifies 
the use of the SAM. Finally we use 

= E» - %(0) - wdTTJdu, (9) 

where Y7 a (ui) is obtained using nonequilibrium PT up to 
second order in u Details of the nonequilibrium RPT were 
published in previous work a 8 i n . The second-order di- 
agram has two sums over Matsubara frequencies. For 
A independent of frequency, one of the sums can be 
done analytically, which simplifies the numerical evalu- 
ation. Explicit expressions for the retarded and lesser 
self-energies are given in the appendix of Ref. 0- It is 
known 21 that dJ7 a /duj = -(3 - n 2 /A)u 2 . 

We begin by analyzing the case u — 1, which corre- 
sponds to the strong coupling (SC) limit U — » 00. As 
in the experiment^ we obtain Tk by a fit of the tem- 
perature dependence of the conductance for V = to 
Eq. © for T/T K < 0.25. In Fig. Q] we compare Eq. 
([3]) with our result. The fit is very good for low T/Tk, 
while for T /Tk — 1, our result lies below the empirical 
curve. Remarkably, this is also the case of the exper- 
imental results i£ This deviation however is outside the 
region of the fit and is irrelevant in the following discus- 
sion. From the fit we obtain VTk = 0.757A. The exact 
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FIG. 1: Conductance shift for V = as a function of tem- 
perature for u — 1. Dashed line corresponds to the empirical 
(EMP) curve Eq. (J3]). 



FIG. 2: (Color online) Scaled conductance shift as a function 
of bias voltage for different temperatures and u = 1. 



results to order T 2 and V 2 can be written in the form 8 

G_ _ ir 2 (l + 2u 2 ) fkT\ 2 
Go 3 \AJ 

_4-3A+(2 + 3^p^ (iQ) 

Note that in the Kondo limit u — ► 1, the coefficients 
are independent of the asymmetry parameter A, in agree- 
ment with recent calculations with the Kondo model.— 

A comparison with the expansion of Eq. ([3]) (for V — 

0) up to second order in T leads to 

kT K 3^(2^ - 1) 
A 7t 2 (1 + 2m 2 )' 1 ' 

which for u = 1 implies T K = 0.746A/fc = 0.949Tjf. The 
small discrepancy with the value kTK — 0.757A obtained 
from the fit is due to the finite temperature interval used 
for fitting. 

Next we calculate the conductance G = dl/dV for 
finite T and V, by numerical differentiation of Eq. ([5]) 
and compare the results with Eqs. and ©. To obtain 
7 we have fitted the current to a polynomial with odd 
powers of V up to F 3 within the range < eV/kTx < 0.4, 
as in the experiments. Inclusion of terms of order V 5 
practically does not modify the results. The resulting 
shift in the conductance (1 — G/Ge(T))/cxv scaled as 
in the experimental work with ay = acr/[l + (7/0; — 

1) c T {T/T K ) 2 } 7 is shown in Fig. H 

From the fit for T — we obtain a = 0.151. This 
agrees with Eq. (JTUJ), which predicts a ratio of the coef- 
ficients of the V and T dependence of the conductance 
[see Eq. (HJ] a = 3/(2vr 2 ) = 0.152 for 5=1, indepen- 
dently of the asymmetry A. . The small discrepancy is 
probably due to numerical errors in the integration near 
the singularities of the integrand^ The resulting values 
of 7 are near 0.75 for 0.13 < T/Tk < 0.4 and increase 



to 0.91 for T/T K = 0.59 and 1.2 for T/T K = 0.79 (not 
shown). These values are larger larger than the experi- 
mentally reported 7 = 0.5 ± 0.1. In addition, while the 
quadratic scaling with V was expected, for most cases 
the observed exponent a = 0.10 ± 0.015 is smaller than 
the value cv|£ M = 3/(2tt 2 ) of the SAM in the SC limit 
U — * 00. Note however that for some of the measured 
systems a approaches this value (Fig. 3 of Ref. for 
Vg ~ —0.195 mV). In addition, in comparison with other 
theoretical predictions for a, 3/ir 2 (Ref. [^), 3/(87r 2 ) (Ref. 
l23h . and 4/7T 2 (Ref. 24), the above value of a lies closer 
to experiment. However, it is clear that the Anderson 
model in the SC limit is not able to reproduce quantita- 
tively the experiment. 

A value of a < agc M can be obtained decreasing u 
[see Eq. (|10p] allowing some degree of intermediate va- 
lence. We explore this possibility, within the SAM, al- 
lowing finite U and therefore u < 1. This means that 
while the average occupation of the dot is kept at the 
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FIG. 3: Conductance shift for V — as a function of temper- 
ature for u — 0.365. Dashed line corresponds to Eq. @. 
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FIG. 4: (Color online) Scaled conductance shift as a function 
of bias voltage for different temperatures and u = 0.365. 

same value {n = 1 in the SAM), some charge fluctu- 
ations with the neighboring configurations is allowed. 
From Eqs. (Q]) and (fTU|) one sees that a = 0.1 implies 
u = 0.365. We have repeated the calculations for this 
value of u. The conductance as a function of temper- 
ature is displayed in Fig. [3] We see that in this case, 
our result lies even closer to the phenomenological curve, 
indicating that the effect of temperature beyond T 2 is 
well described by our approximation (which assumed u 
independent of T and V). Proceeding with the fit as in 
the experiment, we obtain Tk — 1.159A, while Eq. (TTTj) 
gives T K = 1.147A/fc = 1.461T£. 

In Fig. 2] we show the scaled shift in the conductance 
(1 - G/G E (T))/a v [see Eq. ©] for u = 0.365 as a 



function of the applied voltage for several temperatures. 

The corresponding values of 7 are 0.75, 0.57, 0.49, 0.43 
for T/T K =0.17, 0.35, 0.52, and 0.69 respectively. Thus 
except for the smallest temperatures T/Tk < 0.17, in 
this case 7 agrees with the experimental value 7 = 0.5 ± 
0.1. 

The reader might wonder if different physical situa- 
tions in which the Kondo temperature can vary within a 
factor two are consistent with similar values of the renor- 
malized ratio u = U/(irA). In fact, while the Kondo 
temperature decreases exponentially by increasing the ra- 
tio of the bare parameters u — U/ttA, u increases much 
slower,*2ii£ being of course u ~ u for small u (including 
u ~ 0.4) and saturating at u = 1, for u — > 00. 

Finally, we note that if Eq. (1) with arbitrary ex- 
ponents is used to fit the conductance with finite volt- 
age and temperature ranges^ the resulting exponents are 
slightly below 2 in agreement with experiment, 2 

In summary, using nonequilibrium renormalized per- 
turbation theory up to second order in the renormal- 
ized perturbation parameter u for the Anderson model 
in the symmetric case, we have examined the scaling 
behaviour of the conductance, including terms beyond 
those quadratic in temperature and bias voltage. In the 
strong coupling limit, the model predicts an effect of volt- 
age which is 50% higher than observed and the effects of 
terms of order (VT) 2 disagree with experiment. If in- 
stead an important degree of valence fluctuations is al- 
lowed, we obtain good agreement with recent experimen- 
tal results. 
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